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Abstract 


The structure-property relationship is important in understanding molecular 
behaviors and their best-fit areas of applications. 3-(4-hydroxypheny]l) 
prop-2-en-1l-one 4-phenyl Schiff base and some of its derivatives were 
optimized via the density functional theory with Becke three Lee Yang Parr 
correlation and 6-31G* basis set. The molecular properties calculated were 
the energies of the frontier molecular orbitals [highest occupied molecular 
orbital (Exomo), lowest unoccupied molecular orbital (ELumo), energy 
bandgap (Eg), chemical hardness (n), softness (S) and hyperpolarizabilities 
(B)]. The electronic transitions were calculated with the time-dependent 
density functional theory methods, the absorption maxima (Aabs), vertical 
transition energies (AEge), oscillator strengths (f) and molecular orbital 
(MO) components with their percentage contributions were obtained. The 
anti-microbial efficacy of the molecules was tested against Staphylococcus 
aureus aminopeptidase S (AmpS) active site to predict the binding affinities. 
ADMEtox parameters of all the molecules were also investigated. Es values 
ranged from 3.13 to 3.95 eV, B values ranged from 1.45 to 5.81x10-°° esu, 
and their binding affinities ranged from -4.57 to -6.12 kcal/mol, all were 
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more than that of standard drug, streptomycin (-4.31 kcal/mol). The number 
of hydrogen bond donors and hydrogen bond acceptors were ranged from | 
to 2 and 3.75 to 5.25, respectively. Variations observed from the calculated 
molecular properties are the result of varying substituent groups. The 
molecules can be used as nonlinear optical (NLO) materials and also 
showed potential for being effective against Staphylococcus aureus. 


Introduction 


To predict the behavior of molecules correctly and 
ascertain which areas they can best be applied, 
understanding the relationship between structure 
and property is very imperative. A lot of work has 
been reported on how structure influences the 
properties of molecules [1, 2]. Organic Schiff bases 
have been investigated and reported for their 
medicinal applications [3, 4]. Some of them have 
been applied as nonlinear optical (NLO) materials 
and served as a protective shield for the human eyes 
and other sensors from intense laser lights/beams 
[5, 6] while some have been used to inhibit the 
corrosion of metals [7, 8]. Energy band gap 
engineering has been employed by _ synthetic 
chemists while designing novel materials with 
better optoelectronic properties than existing ones. 
Kroon and coworkers [9] highlighted the 
contributing factors influencing molecular energy 
bandgap (Eg) in their work and it has been used as a 
predictive guide in designing and arriving at 
materials with optimal applications. The properties 
of molecules have been predicted via computational 
techniques in the past and in recent times. Density 
functional theory (DFT) has been mostly used in 
predicting the properties of m-conjugated organic 
molecules [10-13]. Properties like polarizability 
and hyperpolarizability [14], geometric parameters 
like bond lengths, bond angles and dihedral angles 
[11], frontier molecular orbitals (FMOs) energies 
[15], global reactivity descriptors [16, 17] amongst 
other properties have been predicted using DFT 
methods. Molecular docking has also been used to 
predict the binding affinities of molecules and their 
mode of action on the receptors of interest [18, 19]. 

3-(4-hydroxyphenyl) prop-2-en-1-one 4-phenyl 
Schiff base and some of its derivatives (Fig. 1) were 
synthesized and reported for their antimicrobial 
potentials [20]. However, the mode of action of 
these molecules on the pathogens screened was not 
ascertained. This study gives an insight into the 
mode of action or mechanism of its derivatives 
(ligands) with an enzyme’s active site. Also, owing 
to the nature of the molecules (extensive 
conjugation), their electronic and NLO properties 


were investigated for probable applications in 
optoelectronic and optical limiting devices. This 
work also predicts the inhibitory action of the 
molecules against Staphylococcus aureus via 
molecular docking of the molecules (ligands) with 
S. aureus aminopeptidase S (AmpS) active sites as 
described by Odintsov and coworkers [21]. It 
further investigates optoelectronic and NLO 
abilities and substituent effects on the molecular 
properties of the molecules. 


N OH 
44 


O=Q) 


R 


Fig. 1 The structure of the molecule under investigation 
where R = OH (f) [21]. In this work, the derivatives were 
uses with substitution at R [(a) R = -H, (b) R = 4-OCHs, (c) 
R = -2-OCH3, 4-OCHs, (d) R = 4-NOz, (e) R = 4-Cl, (f) R= 
4-OH, (g) R = 3-Cl, 4-Cl. 


Methods 


Geometry optimization 


Quantum mechanics calculations were performed 
on the modeled molecules, first by conformational 
search via the molecular mechanics force field 
(MMFF) method to get the true representatives of 
the molecules (lowest energy conformer). After that 
they were subjected to optimization using the DFT 
method with hybrid Becke Three Lee Yang Parr 
(B3LYP) correlation [22] and a polar basis set, 6- 
31G (d) by a restricted hybrid Hartree Fock- DFT 
(HF-DFT) self-consistent field (SCF) calculation 
using Pulay DIIS [23] plus geometric direct 
minimization, all on Spartan 14 computational 
chemistry software [24]. Frequency calculations 
were performed to generate the true ground states 
of the molecules. This level of theory was chosen 
because it has proven to agree with experimental 
trends and findings [14, 17]. The calculated 


molecular properties included the FMOs energies 
(energies of the highest occupied molecular orbital, 
Exomo and that of the lowest unoccupied molecular 
orbital, Etumo), Eg (Eq. 1), hyperpolarizability, 6 
that 1s a measure of the second harmonic generation 
(SHG) efficiency and a third-rank tensor described 
by a3 x 3 x 3 matrix but became 10 components 
instead of 27 components due to Kleimann’s 
symmetry (Eq. 2) [25], hardness, n (Eq. 3) and 
softness, S (Eq. 4) [26]. All these are presented in 
Table 1. 


Eg = ELumo — Exomo Eq. | 


Bere (x 10°? esu) = [(Bxxx + Bxyy + Bxzz)* + (Byyy + Byzz 


ae Does) + (Brzz F Baer ee al Me Eq. 2 
E 

i= = Eq. 3 
1 


Molecular docking 
Protein preparation 


The 3D crystal structure of the human galectin-1 
(PDB ID: 1ZJC) was obtained from the protein data 
bank [27]. AmpS has an active site with two metal 
ions and the protein, 1ZJC was preferred because of 
its high resolution of 1.8 A [21]. A protein 
preparation wizard was used to prepare the 
downloaded protein by assigning charges and bond 
orders and deleting water molecules [28]. OPLS3 
was used for energy minimization [29]. 


Ligand preparation 


The optimized Schiff base molecules were prepared 
using Ligprep in Schrédinger Suite 2017-1 with an 
OPLS3 force field [29]. Epik2.2 in Schr6édinger 
Suite at pH 7.0+2.0 was used to generate the 
ionization states [30]. The prepared ligands were 
docked at the active site of 1ZJC after generating 
receptor grids using OLPS3 force field. XP Ligand 
docking was performed using Glide of Schrédinger 
Maestro v 2017-1 [31]. 
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ADMETYTox screening 


The lead compounds were further screened for 
absorption, distribution, metabolism, excretion and 
toxicity properties using the QikProp program 
(https://www.schrodinger.com/qikprop) and 
Spartan 14 package. These properties are essential 
in testing drugs for their drug-likeness ability and 
save cost involved in bioassay studies. Some basic 
molecular descriptors like molecular weight (Mw), 
number of hydrogen bond donors (donorHB) and 
number of hydrogen bond acceptors (acceptHB), 
etc. were used to ascertain the bioavailability of 
drugs using Lipinski’s rule of five (ROS5) [32, 33]. 
The ROS states that molecules with My<500 have 
good membrane permeability, also the molecules 
must possess donorHB <5 and acceptHB <10 [32, 
34, 35]. Other properties considered were the total 
solvent accessible surface area (SASA), which must 
be between 300.0-1000.0, octanol/water partition 
coefficient (QPlog Po/w), which must be between - 
2.0-6.5, value for serum protein binding 
(QPlogKusa), which must be between -1.5-1.5 and 
other descriptors, together with the binding 
affinities of the molecules (Table 3). 


Results and Discussion 


Electronic properties 


From Table 1, the variation in the Eg values is due 
to different substituent groups attached as seen in 
Fig. Sl. Unsubstituted (R=H), R= 4-OCH3 and R= 
4-OH are almost of similar values while other 
substituted derivatives are of lower values, with R= 
4-NO> being the least. The intramolecular charge 
transfer (ICT) order is as follows: R = 4-NQ2 > R= 
3-Cl, 4-Cl > R= 2-OCH3, 4- OCH3 > R= 4-Cl > R= 
4-OCH3 > R=-H > R= 4-OH. This is because -NO2 
is the strongest electron-withdrawing group (EWG) 
among the substituents followed by —Cl and -OCH3 
while OH is an electron-donating group (EDG). A 
similar trend is observed in their 7 values since 
chemical hardness 1s a direct consequence of energy 
bandgap (Eq. 3), while the S values follow a reverse 


Table 1 Molecular properties of 3-(4-hydroxyphenyl) prop-2-en-1-one 4-phenyl Schiff base and some of its derivatives. 


Molecules Exumo (eV) Exomo (eV) 
R= -H -1.97 -5.90 
R= 4-OCH3 -1.82 -5.74 
R= 2-OCH3, 4- OCH3 -1.72 -5.52 
R= 4-NO2 -2.91 -6.04 
R= 4-Cl -2.10 -5.94 
R= 4-OH -1.85 -5.80 
R= 3-Cl, 4-Cl -2.25 -5.98 


Eg (eV) n (eV) S (eV!) B (x10° esu) 
3.93 1.97 0.509 1.71 
3.92 1.96 0.510 1.59 
3.80 1.90 0.526 2.09 
3.13 hos 0.639 5.81 
3.84 1.92 0.521 4.30 
3.95 1.98 0.506 1.45 
3.73 1.87 0.536 4.50 


Calculated 





260 320 380 440 
UV/Vis spectrum 


Fig. 2 The UV/Vis spectra of the molecules under 
investigation. 


trend (Eq. 4). Since the two FMOs must overlap in 
the middle region of the electron density for 
efficient charge transfer to occur, the ability of the 
molecules to form polar/multipolar components 1s 
influenced by the E, values, with 4-NOz2 most likely 
to possess the highest ability as reported by 
previous researchers [2, 10]. Hence, there is a need 
to investigate the variations of E, values with B 
values. 


Electronic transitions 


NLO-active molecules display good absorption and 
emission properties, the TDDFT calculations were 
used to obtain the electronic vertical singlet 
excitation energies. Their absorption maxima (Aabs), 
oscillator strengths (f), excitation energies 
otherwise known as the vertical transition energies 
(AEce) and percentage contribution of various 
configurations to the excitations were obtained and 
summarized in Table 2 and Fig. 2. These transitions 
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responsible for these absorptions are the m-2* and 
n-z* transitions. It is noteworthy that R = -H, R = 
4-OCH3 and R = 4-OH absorbed around the same 
region while others (2-OCH3, 4-OCHs3, 4-NOz, 4-Cl 
and 3-Cl,4-Cl) bathochromically shifted the Aabs 
significantly, which is attributed to their lower Eg, 
values. These variations and shifts in the absorption 
maxima are indicative that they are correlated with 
the E, values. 4-NO2 has the lowest E, value and 
consequently the highest absorption length, it being 
the best NLO candidate of the studied molecules 
will be confirmed by its / value regarding other 
analogs. 


Nonlinear optical properties 


Molecules exhibiting NLO responses display 
asymmetric polarization brought about by EDGs 
and EWGs once the system is a-conjugated. As it 
has been shown earlier that these groups also altered 
the Es values of the molecules, it is, therefore, 
pertinent to look at the £ values and the relationship 
with Es values. It is evident that the —NO2 
substituted analog has the largest NLO response. 
The order of NLO response is as follows: R = 4- 
NOz > R= 3-Cl, 4-Cl > R= 4-Cl > R= 2-OCHs:, 4- 
OCH; > R = -H > R= 4-OCH3 > R= 4-OH. The NLO 
responses derived show clearly that the greater the 
number and the stronger the EWGs, the greater the 
NLO response. This is because there will be a 
greater push-pull effect in the molecules as EWGs 
are added, which consequently results in larger 
hyperpolarizability values. This correlates with the 
energy band gap earlier discussed that £ values 
increase with the decrease in Eg, values (Fig. 3). 
Therefore, 2-OCH3, 4-OCH3, 4-NO>, 4-Cl and 3- 
Cl,4-Cl derivatives are all better and are identified 
as better choices as NLO materials. However, all 
the studied molecules have higher / values than 


mEg(eV) ™®£ (x10-30 esu) 





R= 4-Cl 


R= 4-OH R= 3-Cl, 4-Cl 


Fig. 3 The relationship/variation of energy bandgap (Eg) with hyperpolarizability (f). 


Table 2 The absorption maxima (Aabs), vertical transition energies (AEge), oscillator strengths (OS) and molecular orbital (MO) 
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components of 3-(4-hydroxyphenyl) prop-2-en-1-one 4-phenyl Schiff base and some of its derivatives. 


Molecules Aabs (nM) AEeg (eV) 
R= -H 349.37 3.55 
291.70 4.25 
R= 4-OCH3 349.96 3.54 
288.54 4.30 
R= 2-OCH3, 4- OCH3 359.26 3.45 
325.30 3.81 
R= 4-NO2 400.64 3.09 
322.28 3.85 
R= 4-Cl 354.02 3.50 
293.00 4.23 
R= 4-OH 348.24 3.56 
286.50 4.33 
R= 3-Cl, 4-Cl 341.93 3.63 
311.82 3.98 


OS (f) 
0.6632 


0.2300 


0.9423 
0.2630 


0.7376 


0.2030 
0.3508 
0.1453 
0.6002 
0.1607 
0.9173 
0.2853 
0.5256 


0.1996 


MO components and contributions 


HOMO -> LUMO 57% 
HOMO-1 -> LUMO 35% 
HOMO-3-LUMO 28% 
HOMO-2-LUMO+1 21% 
HOMO -> LUMO 74% 
HOMO-3-LUMO 46% 
HOMO-LUMO+1 17% 
HOMO -> LUMO 61% 
HOMO-3 -> LUMO 17% 
HOMO-1 -> LUMO 12% 
HOMO-LUMO+1 51% 
HOMO-2-LUMO 20% 
HOMO-1 -> LUMO 85% 
HOMO-3-LUMO 75% 
HOMO -> LUMO 73% 
HOMO-1 -> LUMO 22% 
HOMO-2-LUMO+1 32% 
HOMO-LUMO+1 19% 
HOMO -> LUMO 77% 
HOMO-4-LUMO 71% 
HOMO-1 -> LUMO 47% 
HOMO-?2 -> LUMO 11% 
HOMO-1-LUMO+1 63% 


urea (0.65 x10°°° esu), a standard for organic NLO 
materials [14]. 4-NO2 (5.81 x 10°° esu) has # values 
8.94 times higher than that of urea. All these 
molecules can be used as NLO devices, with the 
best of them being 4-NOz2, R= 4-Cl, R= 2-OCH3, R= 
3-Cl, 4-Cl and 4- OCH:3. 


Antimicrobial potential/docking score 


Molecular docking gives vital information about 
drug/biomolecular interactions, it also probes into 
the mechanisms of action of drugs with the 
receptors to form stable complexes with improved 
efficacy [36]. The information obtained from the 
docking process can be used to suggest the binding 
energy, free energy and stability of complexes. In 
this study, all the molecules showed lower binding 
energy compared to _ the standard drug 
(streptomycin) as presented in Table 3. From the 
order of chemical interaction, all the molecules 
were able to interact with and capture some of the 
reported amino acid residues (His 381, Asp 383 and 
Try 355) [21] within the 4A of AmpS active site. R= 
-H was docked into AmpS active site with a docking 
score of -4.72 kcal/mol as presented in Table 3, and 
the phenol group of R = -H established hydrogen 
bond communication with the reported aromatic 
amino acid Try 355. Arg 31 displayed pi-cation 


HOMO-2-LUMO 17% 


bond to R = -H, which contributes to its stability 
within the AmpS active site as presented in the 
ligand interaction diagram (Fig. S1). R= 4-OCH3 
showed the same order of chemical interaction with 
R = -H but bind higher with -4.57 kcal/mol as 
presented in Fig. S2 with pi-cation interactions 
between the phenyl ring and Arg 31, pi-pi stacking 
with His 381 and hydrogen bonding between the 
hydroxyl group and Tyr 355. R = 2-OCHs3 and 4- 
OCH3 showed the same bond communications with 
above-discussed molecules but for a pi-pi stacking 
with reported aromatic amino acid residue His 381 
with the docking score of -4.59 kcal/mol (Fig. S3). 
The a-x stacking interaction, an attractive non- 
covalent interaction, refers specifically to 
interactions involving aromatic groups containing 7 
bonds. R = 4-NO2 showed no hydrogen bond 
communication with any of the resident amino acid 
residues, but the nitrobenzene group of R = 4-NO2 
communicated with the Glu 106 with the strongest 
non-covalent salt bridge for stabilization in the 
active site, which was established with the docking 
score of -5.47 kcal/mol as presented in Fig. S4. R = 
4-Cl showed no communication with the resident 
amino acid residues but was seen within the 4A of 
the enzyme active site with the lowest docking 
score of -6.12 kcal /mol (Fig. S5). R= 4-OH 
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Table 3 ADME and drug-likeness properties of 3-(4-hydroxyphenyl) prop-2-en-1l-one 4-phenyl Schiff base and some of its 
derivatives against Staphylococcus aureus aminopeptidase S (AmpS) active site. 


Binding affinity Mw 

Molecules (kcal/mol) SASA a 
R= -H -4,72 662.91 327.38 
R= 4-OCH3 -4.57 698.23 357.41 
R= 2-OCH3, 4- OCH3 -4.59 T1dd 387.44 
R= 4-NO2 -5.47 703.18 372.38 
R= 4-Cl -6.12 686.89 361.83 
R= 4-OH -4.87 675.31 343.38 
R= 3-Cl, 4-Cl -4.81 707.43 396.27 
oa 431 804.54 581.57 
(Streptomycin) 


Donor Accpt QPlog QPlog QPlog QPlog 
HB HB Po/w HERG BB Khsa 

1 3.15 4.67 -7.47 -1.06 0.61 

1 4.50 4.76 -7.31 -1.14 0.62 

1 5.25 4.89 -7.24 -1.23 0.64 

1 4.75 3.97 -7.41 -2.23 0.57 

1 a.15 5.16 -7.37 -0.91 0.73 

2 4.50 3.89 -7.34 -1.72 0.41 

1 oe 5.59 -7.28 -0.78 0.84 
16 29.25 -5.93 -5.39 -5.02 -1.90 


Where predicted ICs) value for blockage of HERG K* channels, QPlogHERG = concern below -5; QPlogBB (predicted blood/brain partition 
coefficient) = -3.0-1.2; HB = hemoglobin; SASA = solvent accessible surface area. 


derivative was able to communicate with Tyr 355, 
Glu 106 and Arg 90 by hydrogen bond and was 
stabilized by a pi-cation bond with Arg 31 with the 
docking score of -4.87 kcal/mol (Fig. S6). R = 3-Cl, 
4-C] molecule phenol group communicated with the 
reported aromatic amino acid residue with hydrogen 
bond and stabilized by pi-cation with Arg 31 as 
presented in Fig. S/7/._ Finally, binding 
communication of streptomycin (Fig. S8), a well- 
studied antimicrobial drug for bacteria-derived 
infection. Streptomycin amide and hydroxyl group 
donated hydrogen bond to Glu 106, and His 381 and 
its oxygen group accepted one hydrogen bond from 
Thr 339 for stabilization at the active site but 
showed the highest binding energy of -4.31 
kcal/mol, thus showed it less affinity for the 
Staphylococcus aureus aminopeptidase S (Amp5S). It 
is evident from the study that all the Schiff base 
molecules exert stronger binding affinities for 
inhibition of AmpS as calculated by the glide 
scoring function compared to streptomycin. This 
binding affinity may speak of the anti-microbial 
efficacy of all the molecules. 


ADME analysis 


The ADME (absorption, distribution, metabolism, 
and excretion) properties of all the molecules and 
streptomycin were evaluated to explain their 
pharmacokinetic properties and to determine the 
fitness of molecules as drugs (Table 3). From the 
results, 1t was observed that the blood or brain 
barrier permeability of the standard compound was 
out of the acceptable range, which is very important 
for a drug to pass through those barriers. All the 
molecules under investigation showed QPlogBB 
value range -2.23-1.00, which is better than 
streptomycin (-5.020) where the acceptable range is 
-3.0 to 1.2. The number of hydrogen bond donor 


and acceptor and the solvent-accessible surface area 
(SASA) are acceptable. The predicted ICso value for 
blocking HERG K* channel was very close to the 
acceptable range for streptomycin (-5.390) and 
those of the molecules were within the normal 
range. The predicted octanol or water partition 
coefficient for streptomycin and _ investigated 
molecules were also analyzed. The molecules (-2.0- 
6.5) showed better result than streptomycin (-5.93), 
where the acceptable range is -2.0 to 6.5. The 
predicted serum protein binding values of the 
molecules were within the normal range with the 
value of 0.57-0.84, but streptomycin was out of the 
normal range with a value of -1.9. The normal range 
is -1.5-1.5, thus shows that all the molecules are 
more likely to reach and bind the AmpS active site 
better than streptomycin to elicit the desired 
pharmacologic response. It is evident from the 
ADMETox parameters that all molecules under 
investigation showed less toxic effect compared to 
the streptomycin. 


Conclusions 


The electronic and NLO properties of 3-(4- 
hydroxyphenyl) prop-2-en-l-one 4-phenyl Schiff 
base and its derivatives were predicted using the 
DFT/B3LYP/6-31G*. It was observed that their 
reactivities and NLO responses were enhanced with 
different substituent groups. Their NLO responses 
were compared to that of urea, a standard NLO 
organic material. The binding affinities and 
ADMEtox parameters of the molecules were also 
investigated and compared with those of a standard 
drug (streptomycin). The results show that all the 
molecules will reach and bind the AmpS active site 
better than streptomycin to elicit the desired 
pharmacologic response. It is evident from the 
ADMETox parameters that all molecules under 


investigation showed fewer toxic effects compared 
to streptomycin. Though our observations here are 
still subject to confirmatory procedures for pre- 
clinical and clinical investigations, which best 
probe the atomistic and pharmacology status of a 
drug candidate. 
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